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ABSTRACT 


This  report  considers  the  numerical  determination  of 
the  current  and  charge  induced  on  s tepped-radius  wires  by 
electromagnetic  sources.  A  numerical  technique  is  pre¬ 
sented  which  provides  an  accurate  solution  for  the  current 
and  charge  on  collinear  s tepped-radius  wires.  With  this 
technique  the  detailed  behavior  of  the  current  and  charge 
in  the  vicinity  of  the  step  change  in  radius  is  also 
determined.  A  numerical  technique  employing  the  thin-wire 
theory  approximations  is  also  presented.  This  thin-wire 
approach  requires  no  condition  to  be  enforced  on  the  charge 
at  junctions  of  wires  of  dissimilar  radius  and  yields 
accurate  current  and  charge  distributions  away  from  the 
junction  region.  The  concept  of  a  "jump  condition"  for  the 
charge  on  stepped-radius  wires  is  discussed,  and  a  numerical 
charge  jump  condition  is  computed  and  compared  with  the 
analytical  expression  of  Wu  and  King.  The  use  of  a  charge 
jump  condition  in  other  thin-wire  theory  approaches  which 
require  such  a  condition  is  discussed.  In  particular,  the 
accuracy  of  the  Numerical  Electromagnetics  Code  (NEC), 
which  explicitly  enforces  a  charge  jump  condition  at  junc¬ 
tions  of  wires  of  dissimilar  radii,  is  investigated.  The 
effects  on  numerical  current  and  charge  distributions  of 
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SECTION  I 


INTRODUCTION 

In  recent  years  considerable  effort  has  been  expended  de¬ 
veloping  general  purpose  computer  codes  capable  of  modeling 
complicated  wire  antenna  structures  via  the  method  of 
moments.  These  codes  have  also  been  used,  however,  to  model 
complex  surfaces  by  means  of  grids  of  intersecting  wires  [  1  ]  . 
While  the  use  of  a  wire  grid  to  model  surfaces  is  not  entirely 
satisfactory,  the  simplicity  and  generality  of  the  approach 
have  led  to  its  increasing  usefulness  as  a  tool  for  obtaining 
radar  cross-sections  and  for  assessing  the  performance  of 
antennas  attached  to  or  located  near  such  surfaces. 

The  power  and  flexibility  of  general  purpose  wire  codes 
is  largely  due  to  the  simplicity  of  wire  problems,  which  in 
turn  are  simplified  by  the  use  of  the  so-called  thin-wire 
approximations,  one  of  which  is  that  the  current  flowing  on 
the  surface  of  a  wire  is  assumed  to  be  circumferentially 
invariant  on  each  wire.  However,  when  this  approximation 
is  made,  certain  questions  arise  in  the  formulation  as  to  the 
proper  treatment  of  wire  junctions,  including  junctions  of 
wires  of  dissimilar  radii.  The  usual  procedures  employed  for 
treating  thin  wires  require  us  to  enforce  certain  conditions 
on  the  current  and/or  charge  at  the  junction.  From  the 
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continuity  equation,  one  may  show  that  Kirchoff 's  current 
law  must  be  satisfied  if  it  is  assumed  that  no  charge  can 
accumulate  at  the  junction.  On  the  other  hand,  widely  diver¬ 
gent  opinions  have  existed  as  to  what  conditions  the  charge 
should  satisfy  at  a  wire  junction.  Indeed,  several  different 
charge  junction  conditions  have  been  proposed  [2—6 1  and 
considerable  controversy  has  arisen  as  to  which  is  correct 
[7-10].  In  this  report  we  attempt  to  provide  somewhat  more 
definitive  answers  concerning  the  need  for  and  the  validity 
of  charge  jump  conditions  in  thin  wire  codes. 

In  pursuit  of  the  above  goal,  we  are  led  to  consider 
what  might  be  called  the  canonical  problem  of  two  collinear 
vires  with  dissimilar  radii  which  are  joined  to  form  a 
"stepped-radius"  wire.  Such  a  model  might  apply  to  a  whip 
antenna  mounted  on  a  kingpost  or  other  electrically  thin 
cylinder  which  is  collinear  with  the  antenna  but  of  a  larger 
radius.  To  determine  the  current  distribution  on  such  a 
stepped-radius  structure,  it  is  necessary  to  include  in  the 
model  the  electrical  influence  of  the  step.  However,  it  is 
neither  practical  nor  desirable  in  most  problems  to  actually 
determine  the  detailed  nature  of  the  current  in  the  elec¬ 
trically  small  region  around  such  a  junction.  Instead  one 
would  prefer  to  take  an  approach  similar  to  that  commonly 
used  to  treat  transmission  line  discontinuities.  In  these 


problems,  the  transmission  line  voltages  or  currents  are 


first  determined  at  some  distance  from  the  discontinuity 
and  are  then  effectively  extrapolated  to  the  discontinuity. 
The  discontinuity  in  the  extrapolated  voltage  and/or  current 
is  then  modeled  by  lumped  circuit  elements  which  enable  one 
to  ignore  the  detailed  behavior  of  the  fields  near  the  line 
discontinuity  (i.e.,  the  fields  which  contain  higher  order 
non- t ransmi ss ion  line  modes),  while  preserving  the  electrical 
influence  it  has  on  remote  portions  of  the  structure.  In 
this  study,  we  have  had  a  similar  obj ective--ei ther  the 
development  of  new,  or  the  verification  of  previously  de¬ 
veloped  "lumped"  boundary  conditions  on  the  extrapolated 
values  of  the  current  and/or  charge  at  the  radius  discon¬ 
tinuity  in  the  wire.  The  incorporation  of  such  conditions 
into  thin-wire  codes  should  allow  one  to  accurately  compute 
the  current  and  charge  distribution  away  from  the  radius 
discontinuity  without  detailed  knowledge  of  the  charge  and 
current  in  the  junction  region. 

In  order  to  determine  the  appropriate  boundary  condi¬ 
tions  on  the  extrapolated  current  and  charge  at  the  radius 
discontinuity,  we  first  accurately  determine  the  current  and 
charge  distribution  away  from  the  junction  region.  To  this 
end,  we  present  in  Section  II  a  numerical  technique  which 
accurately  determines  the  current  and  charge  distribution  on 
a  collinear  s t epp ed- r ad i u s  antenna,  including  a  determina¬ 
tion  of  the  detailed  behavior  of  the  current  and  charge  in  the 
electrically  small  region  around  the  junction.  The 
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numerically  determined  current  and  charge  distributions  are 
then  employed  to  obtain  numerically  extrapolated  values  of 
the  current  and  charge  at  the  junction,  from  which  a  charge 
"jump  condition"  may  be  computed  and  compared  with  the 
various  derived  conditions  [2,3]  -  [6,7].  In  Section  II  we 
also  present  numerical  results  based  on  a  thin-wire  technique 
which  does  not  explicitly  employ  a  charge  jump  condition  but 
which  nevertheless  correctly  predicts  the  current  and  charge 
distribution  away  from  the  junction  region. 

In  Section  III  we  compare  results  obtained  via  our  accu¬ 
rate  numerical  model  of  the  s t epp ed- r ad ius  antenna  with 
results  obtained  using  the  Numerical  Electromagnetics  Code 
(NEC)  [11]  which  uses  a  very  "smooth"  basis  set  for  the 
current  representation  and  which  enforces  the  charge  jump 
condition  of  Wu  and  King  [7]  at  discontinuities  in  wire 
radii.  The  accuracy  of  the  code  when  structures  involving 
junctions  of  wires  of  dissimilar  radii  are  modeled  is  dis¬ 
cussed.  Conclusions  and  recommendations  based  on  this  study 
are  presented  in  Section  IV. 
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SECTION  II 


NUMERICAL  INVESTIGATION  OF  THE 
STEPPED- RADIUS  ANTENNA 

In  this  section  we  examine  with  some  care  the  problem  of 
a  s tepped-radius  antenna.  A  numerical  technique  is  presented 
by  which  one  may  accurately  determine  the  current  and  charge 
distributions,  including  the  detailed  behavior  in  the  junc¬ 
tion  region,  for  the  s t epped-radius  antenna.  A  body  of 
revolution  model  is  used  to  represent  the  antenna  and  the 
method  of  moments  [12]  is  applied  to  solve  for  the  unknown 
current  distribution  on  the  antenna.  Since  the  method  of 
moments  is  generally  well-known,  the  numerical  procedure  is 
described  only  briefly,  with  particular  emphasis  on  points 
which  require  special  consideration  for  the  present  geometry. 
Numerical  results  for  current  and  charge  distributions  are 
presented  for  several  selected  cases.  Edge  conditions  for 
the  current  and  charge  in  the  vicinity  of  the  radius  change 
are  also  discussed  and  compared  with  the  numerical  results. 

The  s tepped-r ad ius  antenna  is  also  treated  by  a  thin- 
wire  theory  approach  in  this  section  and  results  are  compared 
with  those  obtained  using  the  body  of  revolution  model. 
Numerical  results  obtained  by  both  approaches  are  used  to 
calculate  numerical  "jump  conditions"  for  the  charge  at  the 


radius  discontinuity  and  these  conditions  are  found  to  he  in 


excellent  agreement  with  the  analytical  jump  condition  of 
Wu  and  King  [7J.  Considerations  in  applying  the  charge  jump 
condition  and  the  appropriateness  of  the  condition  for  junc¬ 
tions  at  bends  and  multiple  wire  junctions  are  also  discussed. 

2 . 1  Analysis  of  the  S tepped-Radius  Antenna  as  a  Body  of 

Revo  lui:ion 

Consider  the  s t ep p e d- r ad ius  antenna  of  Fig.  2.1  which  is 

driven  by  a  delta-gap  voltage  source.  The  dimensions  of  the 

antenna  are  specified  by  the  three  lengths  L^,  L-,  and  L^, 

and  by  the  two  radii  a  ^  and  a  The  variable  s  in  Fig.  2.1 

denotes  arc  length  along  the  surface  of  the  antenna.  At 

this  point  we  stress  the  fact  that  for  many  of  the  cases 

considered  in  this  report  the  scale  of  Fig.  2.1  is  extremely 

deceptive.  For  example,  if  the  radii  shown  in  the  figure 

were  to  remain  unchanged  and  the  total  length  of  a  typical 

antenna  (L^  +  h ^  +  L^)  were  drawn  to  scale,  the  resultant 

antenna  might  well  be  several  hundred  meters  long  for  thin 

antennas  such  that  a.  <  a  <<  L  +  L„  +  L _ .  Thus,  in  order 

1  ^  I  l  5 

to  represent  the  current  and  charge  accurately  in  the  junction 
region  (i.e.,  to  have  several  unknowns  in  the  junction 
annulus  region  and  within  several  radii  of  the  junction)  and 
still  maintain  a  manageable  number  of  unknowns  for  computer 
processing,  the  numerical  technique  must  be  capable  of  treat¬ 
ing  an  enormous  range  of  subdomain  sizes.  The  first  step 
in  this  investigation  is  therefore  to  briefly  describe  the 


6 


numerical  technique  and  verify  that  the  computer  code  can 
indeed  accommodate  the  required  wide  dynamic  range  of  sub- 
domain  sizes. 

In  our  analysis,  the  stepped-radius  antenna  is  driven 
by  a  delta-gap  voltage  source  which  does  not  vary  around 
the  circumference  of  the  wire.  Hence,  the  surface  current 
J(s)  and  the  surface  charge  density  q^Cs)  are  functions  of 
the  arc  length  variable  s  only.  Furthermore,  these  quan¬ 
tities  may  be  converted  to  a  "total"  current  !(s)  and  a 
linear  charge  density  q  ^  (  s  )  simply  by  multiplication  by 
2  it  p  ( s )  ,  where  p(s)  is  the  radius  of  the  wire  as  a  function 
of  s.  The  current  distribution  on  the  stepped-radius  antenna 
is  calculated  via  a  method  of  moments  solution  of  the 
i  n t eg r o-d i f f er en t ia 1  equation 


[-jo)A(r)  -  VD(r)]  •  g  =  -Eini'(r)  •  s, 

r  c  S  , 


(2.1) 


where 
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and  where  E  (r)  is  the  forcing  function.  The  coordinate 
vectors  r  and  r'  locate  the  field  and  source  points,  respec¬ 
tively,  with  reference  to  an  arbitrary  coordinate  origin. 
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For  the  numerical  solution  of  (2.1), 


the  unknown  total 


current  I(s)  and  the  linear  charge  density  q^(s)  are  expanded 
in  a  pulse  basis  set.  The  coefficients  of  the  current  expan¬ 
sion  are  treated  as  unknowns  in  the  problem.  The  coefficients 
of  the  charge  expansion  are  related  to  the  coefficients  of 
the  current  expansion  via  a  finite  difference  approximation 
of  the  continuity  equation 

q  =  j.  y  .  j  .  (2.3) 

SO) 

The  current  and  charge  expansions  are  represented  schemati¬ 
cally  in  Fig.  2.2.  Notice  that  the  current  and  charge 
pulses  are  shifted  with  respect  to  one  another  and  that  a 
half  pulse  of  current  (whose  coefficient  is  zero)  is  used  to 
represent  the  current  at  the  end  of  the  antenna  where  the 
current  is  zero. 

When  the  expansions  for  the  current  and  charge  are 
substituted  into  (2.2)  and  the  resulting  equation  (2.1)  is 
"matched"  at  points  in  the  center  of  the  current  subdomains, 
one  obtains  a  set  of  simultaneous  equations  which  may  be 
solved  to  obtain  values  for  the  coefficients  of  the  current 
expansion  set.  The  coefficients  of  the  expansion  for  the 
charge  are  then  computed  by  a  finite  difference  approximation 
of  the  continuity  equation  (2.3).  More  detailed  explanations 
of  this  numerical  procedure  and  its  features  are  available  in 
[13,  14]. 

i 

i 


Implementation  of  the  numerical  techniques  described  in 
113,  14]  in  a  computer  code  is  fairly  straightforward.  How¬ 


ever,  because  the  dimensions  of  the  s t epp ed- r a d 1  us  antenna 
require  an  extremely  large  dynamic  range  of  subdomain  sizes, 
some  special  precautions  must  be  taken.  First,  one  must 
insure  that  the  integrations  indicated  in  (2.2)  are  performed 
accurately,  particularly  when  the  field  point  is  near  the 
source  region.  Secondly,  several  significant  digits  may  be 
lost  in  the  specification  of  the  subdomain  locations  alone 
due  to  the  use  of  extremely  small  subdomain  sizes  such  as 
occur  near  the  junction  region.  Therefore  if  single  precision 
arithmetic  is  used  in  the  computations  (as  in  this  work), 
one  must  take  special  precautions  to  insure  that  no  additional 
accuracy  is  lost  in  specifying  the  numerical  integration 
quadrature  points. 

For  this  work,  an  existing  general  body  of  revolution 
computer  code  (DBR)  [14]  was  modified  to  treat  the  stepped- 
radius  antenna.  This  modified  code,  herein  referred  to  as 
PEC0,  was  specialized  to  handle  the  large  required  dynamic 
range  in  subdomain  sizes,  to  consider  only  a  circumfer¬ 
entially  independent  current,  and  to  calculate  other  special 
cu r r e n t- d ep end en t  quantities  referred  to  later  in  this  report. 
Verification  that  the  code  can  indeed  handle  a  very  large 
range  of  subdomain  sizes  is  obtained  by  considering  a  simple 
linear  antenna  without  the  step-radius  discontinuity.  Figs. 
2.3  and  2.4  show  the  current  and  charge  distributions, 
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respectively,  for  such  an  antenna  with  radius  ka  ^  =  ka^ 

=  0.000125  and  lengths  k  ( l.  ^  +  I.  ^ )  =  6.25  and  k  I.  ^  =  0.6.  An 
extremely  small  wire  radius  was  chosen  because  the  accuracy 
of  the  self  term  evaluation  is  more  critical  for  smaller 
than  for  larger  radii.  For  this  case,  71  evenly  spaced 
subdomains  have  been  used  to  establish  a  well-converged 
base-line  result.  The  results  shown  in  the  figures  have 
been  verified  by  comparison  with  other  well-established 
wire  codes.  Figs.  2.5  and  2.6  show  the  current  and  charge 
d i s t r i b u t i on s ,  r e s p e c t i ve  ly  ,  obtained  when  unevenly  spaced 
subdomains  are  used,  such  that  the  ratio  of  maximum  to  mini¬ 
mum  subdomain  size  on  the  wire  is  200,000  to  1.  The  subdo¬ 
mains  are  placed  such  that  there  are  two  small  subdomains  of 

—  6 

width  kA  .  =  1.0  *  10  located  near  s/A  =  0.6.  Adjacent 

min 

subdomains  are  progressively  increased  in  width  by  a  factor 

of  2.5  until  the  size  reaches  a  maximum  width  of  kA  =  0.2. 

max 

Subdomains  in  the  vicinity  of  the  feed  point  are  maintained 
at  a  constant  width  of  kA  =  0.1.  If  one  were  to  view  the 
charge  distribution  shown  in  Fig.  2.6  on  an  extremely  ex¬ 
panded  scale,  some  very  minor  anomalous  behavior  would  be 
noted,  but  this  is  not  too  surprising  in  view  of  the  fact 
that  the  ratio  indicated  represents  a  dynamic  range  of  6 
significant  digits,  while  the  computer  on  which  the  code 

runs  maintains  only  8  digit  precision.  One  must  also  keep 

*  t 

in  mind  that  the  charge  is  determined  by  a  finite  difference 


operation  on  the  numerical  current  values. 


ka  =  ka  =0.000125 


Current  distribution  on  a  thin  linear  antenna  obtained  with  evenly 
spaced  subdomains. 


ka  =ka  =0.000125 


Figure  2.4.  Charge  distribution  on  a  thin  linear  antenna  obtained  with  evenly 
spaced  subdomains. 


ka  =ka  =0.000125 


Current  distribution  on  a  thin  linear  antenna  obtained  with 
unevenly  spaced  subdomains. 


Other  tests  have  been  run  to  determine  limitations  of 


the  code.  For  example,  in  one  test  case  the  wire  was  first 
divided  into  equally  spaced  segments  and  then  ten  small  sub- 
domains  were  added  at  a  central  point  on  the  wire  such  that 
there  was  an  abrupt  change  in  subdomain  size  of  ratio  1000 
to  1.  The  current  distribution  obtained  for  this  case  was 
found  to  be  virtually  indistinguishable  from  that  shown  in 
Fig.  2.3.  On  the  other  hand,  the  charge  distribution  exhib¬ 
ited  anomalous  behavior,  but  only  over  the  range  of  the  ten 
small  subdomains.  Elsewhere,  the  charge  distribution  was 
essentially  indistinguishable  from  that  of  Fig.  2.4.  It  is 
likely  that  such  anomalous  behavior  is  due  to  the  fact  that 
the  current  changes  so  slowly  over  the  small  region  that  the 
computation  of  the  charge  by  finite  difference  procedures 
involves  small  differences  and  thus  the  resultant  charge  may 
be  significantly  in  error.  In  general,  however,  the  tests 
which  have  been  performed  lead  to  the  conclusion  that  the 
computer  code  PEC0  does  provide  accurate  representations  for 
both  the  current  and  charge  when  the  subdomain  size  varies 
over  an  enormous  dynamic  range,  as  long  as  the  changes  in 
subdomain  size  are  not  too  abrupt.  In  this  study  we  have 
limited  the  ratio  of  the  larger  to  the  smaller  subdomain 
widths  of  adjacent  charge  subdomains  to  be  no  greater  than 
2.5. 

We  now  consider  a  sample  result  for  a  s tepped-radius 
antenna.  Figs.  2.7  and  2.8  illustrate  the  current  and  charge 
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distributions,  respectively,  on  a  stepped-radius  antenna 
with  the  following  dimensions:  k  a  ^  =  0.01,  k  a  ^  =  0.02, 
kL^  =  3.05,  kL^  =  3.2,  kL^  =  0.6.  In  these  figures  the 
current  and  charge  distributions  are  plotted  as  a  function 
of  the  arc  length  variable  s  and  thus  show  also  the  behavior 
of  the  current  and  charge  in  the  annulus  region  of  the 
junction.  What  appears  to  be  a  single  dashed  line  in  each 
of  the  figures  is  actually  two  dashed  lines  which  delineate 
the  extent  of  the  annular  junction  region.  The  annulus 
is  so  small  for  this  case,  however,  that  the  junction  region 
is  not  visible  in  the  plot.  One  observes  that  the  overall 
distributions  are  much  the  same  as  illustrated  in  Figs.  2.3 
and  2.4,  except,  of  course,  for  the  behavior  of  the  charge 
in  the  vicinity  of  the  junction.  One  notes  upon  inspection 
of  Fig.  2.8,  however,  that  the  overall  "level"  of  the  charge 
distribution  on  each  side  of  the  junction  is  different  if 
one  ignores  the  junction  region  itself  and  merely  extrapolates 
the  sinusoidal-looking  charge  into  the  junction  region.  It 
thus  becomes  clear  how  the  concept  of  a  charge  "jump  condi¬ 
tion"  arises.  One  might,  however,  raise  the  question  as  to 
whether  the  computed  behavior  of  the  current  and  charge  in 
the  vicinity  of  the  junction  is  correct,  or  if  it  represents 
some  numerical  anomaly  which  has  affected  the  entire  charge 
distribution.  We  note,  however,  that  the  approximate  varia¬ 
tion  of  the  charge  near  the  junction  can  be  determined  from 
edge  conditions  [15].  Accordingly,  we  propose  the  following 


Current  distribution  on  a  * tepped-radius  antenna  obtained  by 
treating  the  wire  as  a  body  of  revolution. 


kL  =3.05 


approximate  form  for  the  linear  charge  density  in  the  vicin¬ 
ity  of  the  junction: 


q  (s) 


”  Ba2(Sl  -  s2) (s2  -  s)"1/3 
\  Bp(Sl  -  s)  (s  -  s2)~1/3 

_  Bal  (s  -  SjH®!  -  s  2  >  ~  1  7  3 


S  <  s 


s2  <  S  <  s1 


(  2.  4) 


s  ^  <  s 


where  s  is  the  arc  length  along  the  wire  surface,  p  is  the 

cylindrical  coordinate  variable  on  the  annulus,  s2  is  the 

value  of  the  arc  length  variable  at  the  outer  edge  of  the 

annulus  where  p  =  a2>  s^  is  the  value  of  the  arc  length 

variable  at  the  inner  edge  of  the  annulus  where  p  =  a^,  and 

B  is  an  unknown  constant.  In  justifying  (2.4),  we  note  that 

edge  conditions  [15]  require  that  for  s  near  s2>  the  surface 

charge  density  (proportional  to  the  normal  component  of  E) 

must  vary  as  |s  -  s2|  ,  and  for  s  near  s^,  the  charge 

density  must  vary  as  | s  —  s  ^ |  .  The  linear  charge  density 

“1/3 

must  therefore  vary  as  pjs  -  s2 |  and  p|s  -  SjJ,  respec¬ 

tively.  There  is  one  unknown  constant  associated  with  each 
region,  but  two  of  these  may  be  eliminated  by  noting  that  the 
edge  conditions  require  the  charge  density  near  s^  and  s2  to 
vary  in  the  same  manner  on  either  side  of  the  points  s^  and 
s  2 ,  thus  yielding  (2.4).  Furthermore,  the  linear  charge 
density  is  related  to  the  total  current  by 

«<s>  -  ■  (2-5> 


! 


ill')  that  one  may  determine  the  corresponding  variation  of  Lite 
current  by  simple  integration  in  each  region.  Three  additional 
constants  are  introduced  by  the  integration,  but  two  of  these 
may  be  eliminated  by  enforcing  continuity  of  the  current  at 
the  points  and  s The  resulting  expression  for  the  cur¬ 

rent  f o 1  lows : 


Figs.  2.9  and  2.10  show  the  charge  and  current  distributions, 
respectively,  in  the  vicinity  of  the  junction  for  the  case 
considered  in  Figs.  2.7  and  2.8  but  with  an  expanded  scale. 
The  numerically  obtained  values  are  compared  with  the  ex¬ 
pected  variation,  as  given  by  (2.4)  and  (2.6).  To  obtain 
the  comparison  for  the  charge  (Fig.  2.9)  the  expression  (2.4) 
must  be  "matched"  to  the  numerical  data  at  one  point  in  order 
to  determine  the  unknown  constant  B.  For  the  current  (Fig. 
2.10),  the  same  value  of  B  is  used  and  the  numerical  data 
are  matched  at  one  other  point  to  determine  the  unknown 
constant  C  of  (2.6).  While  the  data-fitting  procedure  does 
not  assure  that  the  magnitudes  (i.e.,  the  undetermined  con¬ 
stants)  of  the  computed  values  are  correct,  the  excellent 
agreement  in  the  variation  of  the  charge  and  current 
verifies  that  the  numerical  values  indeed  satisfy  the  edgp 
conditions  and  hence,  in  all  likelihood,  are  correct. 

Some  measured  data  is  also  available  for  two  stepped- 
radius  monopole  structures  fed  at  the  ground  plane  [8]. 

PEC0  is  also  capable  of  treating  the  s t epped- r ad ius  monopole 
as  a  center  driven  dipole  with  two  step  changes  in  radius  and 
thus  comparison  with  the  measured  data  is  possible.  The 
geometry  may  be  visualized  by  referring  to  Fig.  2.1,  from 
which  a  monopole  with  its  image  is  obtained  simply  by  elimi¬ 
nating  segment  3  and  reflecting  segments  1  and  2  about  the 
feed  point.  Both  monopole  structures  have  three  fixed  di¬ 
mensions:  ka^  =  0.05,  ka^  =  0.02,  and  =  X/4.  In  all  of 

the  measurements,  the  length  of  segment  1  was  fixed  at 
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Figure  2.9.  Comparison  of  the  computed  charge  distribution  in  the  vicini 
the  junction  with  the  approximate  distribution  based  on  edge 
dit ions . 


Figure  2.10.  Comparison  of  the  computed  current  distribution  in  the  vicinity 
of  the  junction  with  the  approximate  distribution  based  on  edge 
conditions . 


either  L  j  -  A  /  4  or  L  =  A  /  2  .  Figs.  2.11  and  2.12  nmip.irc  t  lu¬ 
mens  tired  data  with  current  and  charge  distributions,  respec¬ 
tively,  calculated  by  P  E C 0  for  the  case  L  -  A / 4 .  Figs.  2.11 
and  2.14  illustrate  the  comparison  for  the  case  Lj  =  A/2. 

For  both  cases  the  calculated  and  measured  current  distribu¬ 
tions  agree  in  shape.  The  magnitudes,  however,  disagree  by 
a  factor  of  nearly  1.25  but  in  opposite  directions  in  the  two 
cases.  On  the  other  hand,  the  calculated  and  measured  charge 
distributions  are  in  excellent  agreement  for  both  cases.  No 
definitive  explanation  for  the  magnitude  differences  in  the 
calculated  and  measured  current  is  offered  at  this  time. 
However,  there  is  a  strong  possibility  that  the  differences 
arise  due  to  the  difficulty  in  determining  an  accurate  cali¬ 
bration  factor  for  the  measurement  probes  [16].  Further 
data  is  presented  in  [8],  however,  which  demonstrates  agree¬ 
ment  between  calculated  and  measured  data  for  the  case  when 
L  =  A/4.  This  comparison  is  given  in  Fig.  2.15  which  shows 
that  results  obtained  with  PEC0  are  in  excellent  agreement 
with  the  measurements.  However,  the  measured  data  in  Figs. 
2.11  and  2.15  supposedly  apply  to  the  same  case,  so  that  an 
inconsistency  appears  in  the  measured  data.  Furthermore, 
since  the  measured  charge  agrees  with  the  computed  charge 
and  since  the  current  can  be  obtained  from  the  charge  by 
integration,  beginning  from  the  end  of  the  monopole  where  the 
current  is  known  to  vanish,  the  charge  and  current  measure¬ 
ments  shown  in  Figs.  2.11-2.15  are  also  inconsistent.  It  is 
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Comparison  of  calculated  current  distribution  with  measured  data 
for  a  stepped-radius  monopole  of  height  L.+L=A/2. 
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Comparison  of  calculated  current  distribution  with  measured  data  for 
a  stepped-radius  monopole  of  height  L.+L  =  0. 75X . 


possible  Chat  the  measured  data  shown  in  Fig.  2.15  have  been 
scaled  to  demonstrate  the  agreement  in  shape  between  the 
measured  and  the  numerical  data  available  at  that  time  L 1 6  J - 
Unfortunately,  if  this  is  the  case,  mention  of  the  data 
scaling  was  omitted  from  the  text  of  [8].  We  conclude,  how¬ 
ever,  based  on  the  excellent  agreement  between  the  computed 
and  measured  charge  and  on  other  evidence  presented  in  this 
section,  that  the  computer  code  PEC0  does  indeed  provide 
accurate  current  and  charge  distributions  on  stepped-radius 
antenna  structures — including  showing  the  detailed  behavior 
of  the  current  and  charge  in  the  junction  region. 

2 . 2  Thin-Wire  Analysis  of  the  S t epped-Radius  Antenna 

The  analysis  of  the  stepped-radius  antenna  by  means  of 
thin-wire  theory  differs  from  the  body  of  revolution  analysis 
in  two  respects.  First,  the  current  and  charge  on  the  annular 
region  of  the  junction  are  completely  ignored  and  secondly, 
the  so-called  reduced  kernel  is  employed  in  the  integral 
equation.  The  reduced  kernel  is  obtained  by  approximating 


as 


■  j  k [ ( £  -  l')2  +  a2] 


LU  -  V)2  +  a2] 
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where  £  is  defined  in  Fig.  2.1  and  a  is  the  radius  of  the 
wire.  For  the  s tepped-radius  wire,  however,  the  radius  a 
has  two  distinct  values,  and  one  might  ask  whether  the  value 
of  a  to  be  used  is  the  radius  at  the  "match  point"  or  the 
radius  associated  with  the  source  region.  The  answer  depends 
on  the  physical  interpretation  of  (2.8).  One  common  approach 
is  to  interpret  (2.8)  as  effectively  replacing  the  wire  by  a 
filament  of  current  flowing  ctnwn  the  axis  of  the  wire  with 
the  match  point  placed  on  the  wire  surface.  With  this  inter¬ 
pretation,  one  should  obviously  use  the  radius  of  the  wire 
at  the  match  point.  On  the  other  hand,  one  may  also  view  the 
problem  in  an  extended  boundary  condition  sense  [17].  This 
approach  would  imply  that  the  match  point  is  on  the  wire 
axis,  while  the  current  flows  along  the  wire  surface.  The 
appropriate  radius  to  use  in  (2.8)  in  this  case  would 
therefore  be  the  radius  associated  with  the  source  region. 

The  extended  boundary  condition  argument  and  approach  is 
used  in  the  analysis  presented  here  to  avoid  an  ambiguity  in 
the  choice  of  the  match  point  location  at  the  junction. 

As  we  have  said,  in  a  thin-wire  analysis  the  charge  and 
current  on  the  annular  region  are  ignored.  One  may  visualize 
the  placement  of  the  subdomains  in  the  thin-wire  analysis 
then,  simply  by  letting  the  width  of  the  annulus  region  in 
Fig.  2.2  be  equal  to  zero.  Thus  a  single  current  subdomain 
is  split  at  the  junction  while  full  pulses  of  charge  are 
placed  on  each  side  of  the  step.  The  current  subdomain  split 
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hy  the  junction  has  associated  with  it  a  single  unknown  coef¬ 
ficient  while  the  full  pulses  of  charge  residing  on  each  side 
of  the  radius  change  are  completely  independent  of  one  an¬ 
other.  Note  that  the  two  portions  of  the  current  pulse  s  p 1  i  t 
across  the  junction  are  associated  with  different  radii 
according  to  which  section  of  the  wire  they  reside  on.  Also 
note  that  we  need  not  explicitly  enforce  any  condition  on 
the  charge  at  the  junction,  the  only  condition  explicitly 
enforced  at  the  junction  being  that  of  current  continuity. 

We  note,  however,  that  the  testing  procedure,  which  in  this 
approach  requires  that  the  average  tangential  electric  field 
vanish  over  each  current  subdomain,  effectively  (in  the  limit 
of  decreasing  subdomain  size)  enforces  continuity  of  the 
scalar  potential  at  the  junction.  We  emphasize  that  t_h  is 
condition  is,  in  reality,  the  "correct  junction  condition" 
and  hence  we  should  anticipate  reasonable  performance  with 
the  thin-wire  code. 

The  resultant  current  and  charge  distributions  computed 
with  2  7  unknowns  by  this  method  for  the  same  stepped-radius 
antenna  as  that  of  Fig.  2.7  are  shown  in  Figs.  2.16  and  2.17, 
respectively.  The  current  distribution  is  virtually  identical 
to  that  obtained  from  the  body  of  revolution  computer  code 
PEC0  (Fig.  2.7).  Note  that  owing  to  the  sparseness  of  the 
sample  points  in  the  neighborhood  of  the  junction,  the  charge 
distribution  shown  in  Fig.  2.17  does  indeed  appear  to  have  a 
jump  discontinuity  there.  The  computed  values,  however,  are 
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Figure  2.16.  Current  distribution  on  a  stepped-radius  antenna  computed  by 
thin-wire  analysis. 


Imaginary 


virtually  identical  with  the  charge  distribution  computed 
with  the  body  of  revolution  code  (Fig.  2.8). 

Consider  next  a  stepped-radius  antenna  with  the  follow¬ 
ing  dimensions:  k  a  ^  =  0.001,  k  a  ^  =  0.05,  k  L  ^  =  3.05, 
kl^  =  3.2,  and  kL^  =  0.6.  Note  that  the  structure  has  a 
radius  discontinuity  of  ratio  50  to  1.  The  current  and 
charge  distributions  obtained  using  PEC0  are  presented  in 
Figs.  2.18  and  2.19,  respectively.  Figs.  2.20  and  2.21  show 
the  current  and  charge  distributions  obtained  using  the  thin- 
wire  theory  approach  with  26  unknowns.  A  comparison  of  the 
appropriate  figures  shows  that  the  results  for  the  current 
are  in  excellent  agreement.  The  charge  distributions  also 
agree  very  well.  In  Figs.  2.22  and  2.23  the  number  of  un¬ 
knowns  used  in  the  thin-wire  approach  has  been  increased  to 
A3.  The  current  and  charge  distributions  obtained  are  essen¬ 
tially  the  same  as  shown  in  Figs.  2.18  and  2.19.  In  Fig. 
2.23,  however,  one  notes  that  the  charge  in  the  vicinity  of 
the  junction  tends  to  show  the  proper  edge  behavior  (compare 
with  Fig.  2.19)  even  though  the  annulus  region  has  been 
completely  ignored  In  the  thin-wire  approach.  This  behavior 
has  been  verified  for  many  different  wire  radii  and  radius 
jump  ratios,  and  is  further  demonstrated  in  Section  III. 

Thus,  Figures  2.18  through  2.23  demonstrate  that  the  thin- 
wire  approach  presented  in  this  section  can  be  used  to  treat 
wires  with  step  changes  in  radius— even  very  large  changes-- 
without  the  explicit  e n f o r c e men t  of  any  conditions  on  the 
charge  at  the  junction. 
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(PEC0)  on  a  stepped-radius  antenna 
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Figure  2.21.  Charge  distribution  (thin-wire  approach)  on  a  stepped-radius 
antenna  with  a  50  to  1  change  in  radius. 
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Figure  2.22.  Current  distribution  (thin-wire  approach;  increased  number 
unknowns)  on  a  stepped-radius  antenna  with  a  50  to  1  change 
in  radius. 
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Figure  2.23.  Charge  distribution  (thin-wire  approach;  increased  number  of 
unknowns)  on  a  stepped-radlus  antenna  with  a  50  to  1  change 
in  radius. 


2 . 3  Computation  of  the  Chargo  "Jump  Condition" 

In  this  section  numerically  derived  charge  jump  condi¬ 
tions  are  computed  for  various  wire  radii  and  changes  in 
radii  and  are  compared  with  the  analytical  jump  condi ton  ol 
Wu  and  King.  We  compute  the  charge  jump  condition  for  a 
specific  structure  by  ignoring  the  actual  behavior  of  the 
charge  in  the  junction  region  and  instead  extrapolate  to  the 
junction  using  values  of  current  and  charge  computed  on  both 
sides  of  the  junction  at  some  distance  away  from  it.  This 
process  is  analogous  to  that  of  determining  an  equivalent 
circuit  for  a  discontinuity  at  a  terminal  plane  in  a  single 
mode  waveguide — from  measurements  of  standing  waves  away 
from  the  discontinuity,  we  deduce  the  reflection  coefficient 

(and  hence  the  discontinuity  of  the  modal  coefficients) 

a  t  the  terminal  plane,  thus  determining  the  equivalent  cir¬ 
cuit.  This  mode  alone,  however,  would  not  be  useful  for  des¬ 
cribing  the  true  fields  in  the  region  of  the  discontinuity 
because  of  the  presence  of  higher  order  modes  in  the  terminal 
plane  region. 

We  calculate  charge  jump  conditions  using  charge  distri¬ 
butions  obtained  with  both  the  body  of  revolution  code  l’EC0 
and  with  the  thin-wire  approach  described  in  the  previous 
section.  Calculation  of  the  charge  jump  condition  from  the 
thin-wire  approach  can  be  used  to  further  validate  the  thin- 
wire  approach  over  a  wide  range  of  wire  sizes  and  jump  ratios. 

The  calculation  of  the  jump  condition  relies  on  the  fact 
that  away  from  the  junction  region  (and  also  from  wire  ends 
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and  feed  points)  both  the  current  and  the  charge  on  the 
antenna  are  essentially  sinusoidal.  If  the  sinusoidal  varia¬ 
tion  is  independently  determined  on  each  side  of  the  junc¬ 
tion  and  the  resulting  variation  is  extrapolated  to  the 
junction  region,  one  obtains  a  total  current  distribution 
which  has  a  slope  discontinuity  at  the  junction  and  hence  a 
charge  distribution  which  exhibits  a  step  function  jump 
there.  The  extrapolated  charge  variation  in  the  vicinity 
of  the  junction  may  therefore  be  expressed  as 

f~  A^cos(ks)  +  BjSin(ks)  ,  <  s 


q (s)  *  S 


(2.9) 


^A2cos(ks)  +  B^siniks)  ,  s  < 


where  k  is  the  free  space  wavenumber,  s,  s ^ ,  and  s2  are  de¬ 
fined  as  in  Section  2.1,  and  A^ ,  A,,  ,  B^,  and  B2  are  unknown 
constants.  To  determine  the  constants  in  (2.9)  it  is  neces¬ 
sary  to  find  a  region  on  each  side  of  and  away  from  the 
junction,  over  which  the  computed  current  (and  hence  the 
charge)  is  essentially  sinusoidal.  These  regions  are  located 
by  considering  the  expression 


P 


100 


b’*&) 


I(s) 


I  k2I , 


max 


(2.10) 


The  quantity  P  may  be  used  as  a  gauge  of  the  percentage 
deviation  from  sinusoidal.  The  current  is  known  only  as 
. umerical  values,  so  a  finite  difference  approximation  may 
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be  used  to  calculate  the  derivatives  in  (2.10).  I  1 

ina  x 

represents  the  largest  value  of  the  current  magnitude  on  the 
structure,  ignoring  the  region  near  the  feed  point.  The 
points  from  which  the  charge  is  extrapolated  are  found  by 
assuming  a  very  small  value  for  a  quantity  P  and  then 
searching  for  the  points  nearest  the  junction  for  which 
P  <  P  •  The  value  of  P  is  increased  progressively 
until  such  points  can  be  found.  In  practice,  such  regions 
are  generally  found  for  P  <15  (per  cent)  .  The  numerical 

charge  distribution  is  then  equated  to  the  appropriate  repre¬ 
sentation  in  (2.9)  at  two  points  ir.  each  region  to  determine 
the  constants.  The  extrapolated  ch  rge  at  the  junction  is 
then  f  ound  to  be 

9j  =  A^cos(ks^)  +  B^sin(ks^)  (2.11a) 

9  2  =  A^cosCks.,)  +  B^sin  (ks^ )  ,  (2.11b) 

where  q^  and  q2  represent  the  charge  at  the  junction  on  the 
wire  of  smaller  radius  and  larger  radius,  respectively. 

An  example  of  the  results  this  extrapolation  procedure 
is  illustrated  in  Fig.  2.24  for  the  charge  distribution  shown 
in  Fig.  2.19.  The  figure  shows  an  expanded  plot  of  the 
numerically  computed  values  in  the  junction  region.  The 
solid  lines  represent  the  linear  charge  density  (2.9)  after 
the  constants  have  been  computed  by  "matching"  the  data  at 
two  points  on  each  side  of  the  junction.  The  points  at 
which  the  solid  lines  intersect  the  two  dashed  lines 
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ure  2.24.  Charge  distribution  (PEC0)  in  the  vicinity  of  the  junction  on 
stepped- radius  antenna  with  a  50  to  1  change  in  radius.  The 
solid  lines  represent  the  "extrapolated"  charge. 


(delineating  the  annulus  region)  represent  the  real  and 
imaginary  parts  of  q  and  q  .is  defined  by  (  1  .  I  1  )  .  A  number 
of  similar  calcualtions  have  been  made  and  the  resulting  ratio 
q  l  /  q  ?  compared  w  i  t  It  the  formula  of  Wu  and  King  17]  in 

fig.  2.25  for  several  wire  radii  and  various  ratios  of  radii 
at  the  discontinuity.  for  reference,  the  formula  of  Wu  and 
King  is  given  as 


q 

q 


2 


(2.12) 


where  'j  =  0.5772  is  Euler's  constant.  One  observes  from 

Fig.  2.25  that  the  numerical  computations  of  both  the  body 
of  revolution  code  and  the  thin-wire  code  are  in  very  good 
agreement  with  the  formula  of  Wu  and  King.  The  only  excep¬ 
tion  appears  on  the  curve  for  k a  =  0.01,  where  it  is  noted 
that  for  a  ^  /  a  7  =  0.02,  the  point  computed  by  thin-wire 
analysis  is  significantly  in  error.  However,  in  this  case, 
the  radius  of  the  thicker  wire  (  k  a  2  =  0.5)  is  well  beyond 
the  limit  for  which  thin-wire  theory  is  applicable.  We 
comment  that  all  data  were  obtained  for  the  wire  length 
dimensions  k  L  ^  =  3.05,  k  L  2  =  3.2,  and  k  L  ^  =  0.6,  since  this 
results  in  an  approximate  charge  maximum  in  the  vicinity  of 
the  junction,  thus  improving  the  sensitivity  of  the  compu¬ 
tations.  Results  obtained  for  several  other  wire  lengths. 
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Figure  2.25.  Comparison  of  numerically  computed  charge  "jump  condition" 
for  stepped-radius  wires  with  theory  of  Wu  and  King. 
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however,  where  a  charge  maximum  does  not  occur  near  the 
junction,  have  also  confirmed  agreement  of  the  computations 
with  the  formula  of  Wu  and  King.  We  mention  briefly  that  a 
similar  procedure  to  that  above  for  extrapolating  the  charge, 
has  been  used  to  verify  that  the  extrapolated  current  is 
continuous  at  the  junction  in  both  the  body  of  revolution  and 
in  the  thin-wire  theory  approach. 

These  computations  verify  that  the  charge  jump  condition 
of  Wu  and  King  is  correct  for  the  s t ep p ed- r ad i us  structure 
illustrated  in  Fig.  2.1.  Time  and  resources  have  not  permit¬ 
ted  us  to  attempt  to  calculate  charge  jump  conditions  for 
more  complex  dissimilar  radii  wire  structures  such  as  multi¬ 
ple  wire  junctions  or  wires  bent  at  a  radius  discontinuity. 
However,  it  can  easily  be  shown  (see  Appendix)  that  the  Wu- 
King  condition  is  approximately  equivalent  to  requiring 
continuity  of  scalar  potential  at  a  junction.  Hence  estab¬ 
lishing  the  accuracy  of  the  condition  for  the  stepped- radius 
cylinder  lends  a  great  deal  of  support  to  its  applicability 
to  these  more  general  junctions  where,  again,  the  correct 
junction  conditions  are  that  (1)  Kirchoff's  current  law  be 
satisfied  and  that  (2)  scalar  potential  be  continuous  at 
each  junction. 

Because  of  the  difficulty  involved,  it  is  not  likely 
that  a  detailed  investigation  of  bent  wire  junctions  or 
multiple  wire  junctions  at  the  detail  we  have  considered  here 
will  be  carried  out  in  the  near  future.  Nevertheless,  the 
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success  of  our  thin-wire  code  in  modeling  the  stepped-radius 
junction,  leads  us  to  suggest  a  less  ambitious  investigation-- 
that  of  establishing  whether  or  not  the  t h i n- w ire  approach 
described  here  yields  results  satisfying  the  Wu-King  condition 
on  such  structures.  In  applying  the  thin-wire  approach, 
Kirchoff's  current  law  would  be  enforced  at  the  junction,  but 
no  condition  on  the  charge  would  be  required.  Our  conjecture 
is  that  the  proposed  thin-wire  approach  would  yield  charge 
distributions  in  good  agreement  with  the  Wu-King  condition 
because  both  approaches  employ  an  approximate  enforcement  of 
continuity  of  scalar  potential  at  the  junction.  Such  agree¬ 
ment  would,  in  our  opinion,  lend  substantial  credence  to 
the  use  of  both  thin-wire  theory  and  the  Wu-King  condition 
to  handle  such  configurations. 


SECTION  III 


APPLICATION  OF  THE  NUMERICAL  ELECTROMAGNETICS 
COMPUTER  CODE  TO  STEPPED-RADI  'JS  JUNCTION’S 

The  Numerical  Electromagnetic  r.  Code  (NEC)  is  a  general 
purpose  computer  code  capable  of  performing  a  variety  of 
tasks  related  to  electromagnetic  radiation  and  scattering 
problems  111].  NEC  is  based  on  a  method  of  moments  solution 
I  12]  of  integral  equations  for  the  currents  induced  on  a 
structure  (or  structures)  by  electromagnetic  sources  or 
incident  fields.  The  code  is  capable  of  treating  wire  struc¬ 
tures  and/or  smooth  surfaces  (by  wire-grid  modeling  or  the 
Magnetic  Field  Integral  Equation).  In  this  report  we  concern 
ourselves  only  with  a  particular  aspect  of  the  code--the 
modeling  of  junctions  of  wires  of  dissimilar  radii.  This 
single  aspect  is  of  interest  because  NEC  explicitly  employs 
the  charge  jump  condition  (2.12)  of  Wu  and  King  [7]  at  wire 
junctions.  In  Section  II,  we  have  numerically  verified  that 
this  charge  jump  condition  is  indeed  a  "correct"  condition 
to  enforce  at  the  junction  of  two  collinear  wires  and  thus 
we  expect  NEC  to  provide  accurate  results.  On  the  other 
nand,  some  aspects  o f  th in-wire  modeling  indicated  in  the 
previous  section  lead  one  to  expect  that  the  explicit  enforce¬ 
ment  of  (2.12)  might,  under  some  conditions,  result  in 
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erroneous  solutions.  Specifically,  wi1  refer  to  the  fact  that 
ttic  current  and  charge  obtained  via  the  thin-wire  approach 
of  Section  II  tend  ft)  show  the  true  edge  behavior  —  not  the 
extrapolated  behavior — as  the  number  i  f  unknowns  is  increased 
in  the  junction  region. 

In  this  section  we  address  the  question  of  whether  the 
explicit  enforcement  by  NEC  of  the  charge  jump  condition 
(2.12)  at  junctions  of  col  linear  wires  leads  to  accurate 
solutions  for  the  current  and  charge  away  from  the  junc Li  on 
region.  A  very  brief  description  of  the  wire  modeling  tech¬ 
nique  used  in  NEC  is  first  presented  as  background  (for 
further  details,  one  should  see  I.  111).  Numerical  results 
obtained  with  NEC  are  then  compared  with  results  obtained 
using  PEC0  for  various  structures  and  the  data  are  discussed. 

3 . I  Wire  Modeling  in  NEC 

The  NEC  program  models  arbitrary  wire  structures  as  a 
set  of  straight  wire  segments.  for  each  of  these  wire  seg¬ 
ments  one  of  three  thin-wire  approximations  for  the  kernel 
of  the  integral  equation  is  used.  For  most  cases  a  reduced 
kerneL  corresponding  to  (2.8)  is  used,  where  the  current 
is  assumed  to  be  a  filamentary  current  flowing  along  the 
axis  of  each  segment  and  the  match  point  is  assumed  to 
reside  on  the  wire  surface.  For  thick  wires  an  extended 
thin-wire  kernel  fill  can  he  selected  by  the  user.  In  this 
case  the  current  is  assumed  to  flow  on  the  surface  of  the  wire 
and  the  integral  (2.7)  is  approximated  hv  a  two  term  series. 


Finally,  whenever  wire  segments  are  separated  by  large 
distances,  fields  are  calculated  by  treating  the  current 
segment  as  an  infinitesimal  current  element  with  the  ap¬ 
propriate  dipole  moment.  The  current  element  in  this  case 
is  located  at  the  center  of  the  source  segment. 

The  method  of  moments  is  applied  by  choosing  an 
appropriate  set  of  testing  functions  and  basis  functions  for 
the  current  on  the  wire  segments.  Substitution  into  the 
integral  equation  for  the  structure  under  consideration 
results  in  a  set  of  simultaneous  equations  which  can  be 
solved  to  determine  the  coefficients  of  the  current  basis 
functions.  In  the  development  of  the  NEC  program,  the  basis 
functions  used  to  represent  the  current  on  wires  were  chosen 
so  as  to  yield  a  highly  efficient  and  accurate  solution. 

Since  the  current  distribution  on  wire  structures  is  smooth 
(except  at  junctions),  a  very  "smooth"  set  of  basis  functions 
was  chosen  for  NEC.  Furthermore,  because  of  the  accuracy 
with  which  these  smooth  basis  functions  can  represent  a 
smoothly-varying  current,  a  set  of  delta  functions  could  be 
used  for  the  testing  functions,  resulting  in  a  collocation 
or  point-matching  solution  technique.  The  use  of  very 
smooth  basis  functions  generally  results  in  improved  conver¬ 
gence  over  pulse  or  piecewise  linear  current  expansions  and, 
indeed,  we  have  found  the  NEC  program  to  be  a  highly  efficient 
and  accurate  code.  It  is,  however,  the  use  of  smooth  basis 
functions  which  necessitates  the  use  of  an  equivalent  charge 
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jump  condition  to  replace  the  actual  non- smooth  (in  fact, 
singular)  variation  of  the  c  it  urge  at  a  junction. 

The  current  expansion  in  NEC  consists  of  three  terms 
defined  on  each  segment--a  constant,  a  sine,  and  a  cosine 
term.  The  trigonometric  terms  have  the  form  of  standing 
waves  employing  the  free  space  propagation  constant.  The 
total  current  on  each  segment  thus  takes  the  form 

I.U)  =  A.  +  B.  sin  [  k  ( £  -  £.)]  +  C.  cos  [  k  (  £  -  £.)],  (3.1) 

I  £  -  £  .  I  A  .  /  2  , 

1  l  1  t 

where  £  is  defined  in  Fig.  2.1,  £  is  the  value  of  1.  at  the 
center  of  segment  i,  and  A  ^  is  the  length  of  segment  i. 

There  are  three  unknown  constants  (A^,  B^,  C^)  associated 

with  the  current  expansion.  One  of  these  constants  is 
eliminated  by  requiring  that  Kirchoff's  current  law  be  sat¬ 
isfied  at  junctions.  (Between  two  adjacent  wire  segments, 
the  requirement  simply  reduces  to  a  continuity  of  current 
requirement.)  A  second  constant  is  eliminated  by  requiring 
a  condition  on  the  charge  to  be  satisfied  at  junctions.  The 
condition  used  in  NEC  is  the  Wu-King  charge  jump  condition 
given  by  (2.12).  For  two  adjacent  segments  of  the  same1 
radii,  this  condition  corresponds  simply  to  enforcement  of 
continuity  of  charge  at  the  junction.  The  actual  implementa¬ 
tion  of  these  two  conditions  in  NEC  is  somewhat  involved  for 


general  wire  structures  and  the  details  can  be  found  in  [11]. 
The  final  result  of  applying  the  current  and  charge  condi¬ 
tions,  however,  is  that  only  one  unknown  constant  related 
to  the  current  amplitude  on  each  segment  remains  to  be  evalu¬ 
ated  via  the  method  of  moments.  With  this  brief  background, 
we  now  proceed  to  discuss  the  results  obtained  using  NEC. 

3 . 2  Numerical  Results  and  Discussion 

In  this  section  we  wish  to  determine  if  NEC's  explicit 
enforcement  of  the  charge  jump  condition  (2.12)  yields  ac¬ 
curate  current  and  charge  distributions  away  from  junction 
regions.  The  most  reliable  way  to  accomplish  this  appears 
to  be  to  make  direct  graphical  comparisons  of  results 
computed  using  NEC  and  using  the  baseline  code  P  E  C  0  .  Simpler 
comparisons  such  as  of  the  input  idmittance  computed  by 
both  codes,  have  been  found  to  be  less  reliable  because  the 
differences  between  the  codes  in  numerical  modeling  of  the 
source  result  in  slightly  different  admittance  values.  Fol¬ 
lowing  the  presentations  of  the  graphical  comparisons,  an 
attempt  is  made  to  reduce  the  principal  results  of  the  data 
to  a  more  compact  form. 

First,  we  describe  the  format  of  the  data  presented  in 
this  section.  For  all  of  the  figures  illustrating  current 
and  charge  distributions,  results  obtained  using  NEC  are 
plotted  with  different  symbols  for  the  real  and  imaginary 
parts.  Results  obtained  using  the  baseline  code  PEC0  are 
plotted  as  a  solid  line.  In  general,  the  number  of  unknowns 
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used  in  NEC  is  not  specified,  but  this  can  be  easily  deter¬ 
mined  from  the  figures  since  each  symbol  represents  the 
value  of  tiie  current  at  the  center  of  a  subdomain.  Of 
more  importance  is  the  subdomain  width  in  terms  of  wire- 
radii,  which  may  be  different  on  opposite  sides  of  the  junc¬ 
tion.  This  information  is  given  in  most  plots  as  A  .  =  C  a  ^  , 

where  A.  is  the  subdomain  size  and  a.  is  the  radius  of  the 
1  1 

i  wire  (Fig.  2.1)  and  C  is  the  proportionality  factor 
relating  them.  We  note  that  the  subdomain  sizes  indicated 
in  the  figures  refer  only  to  the  two  subdomains  on  either 
side  of  the  junction. 

We  begin  our  comparisons  by  investigating  s t epped- rad i us 
antennas  with  various  radii  and  ratios  of  radii,  but  with 
fixed  lengths.  The  lengths  of  the  structure  considered  are 
(c.f.  Fig.  2.1)  kL  ^  =  3.05,  kL.,  =  3.2,  and  kL^  =  O.b,  chosen 
so  as  to  keep  the  source  a wav  from  the  junction  region  and  to 
result  in  a  charge  maximum  in  the  neighborhood  of  the  junc¬ 
tion.  We  first  consider  several  cases  for  which  the  radius 
of  the  smaller  wire  (section  1)  is  ka  ^  =  0.0001.  Figs.  3.1 
and  3.2  show,  respectively,  the  current  and  charge  distribu¬ 
tions  on  such  an  antenna  with  k a  =  0.000125.  The  radius 
change  in  this  case  is  very  small  and  the  agreement  between 
PEC0  and  NEC  is  excellent.  Figs.  3.3  and  3.4.  show,  respec¬ 
tively,  the  current  and  charge  distributions  with  k  a  ,, 
increased  to  0.001,  representing  a  10  to  1  change  in  radius. 


We  note  that  there  is  a  small  difference  in  the  current  on 
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Comparison  of  calculated  current  distribution  (NEC)  with  baseline 
code  results  for  a  very  thin  antenna  having  a  radius  change  of  ra 
1.25  to  1.0. 


Comparison  of  calculated  charge  distribution  (NEC)  with  baseline 
code  results  for  a  very  thin  ^..tenna  having  a  radius  change  of 
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Figure  3.3.  Comparison  of  calculated  current  distribution  (NEC)  with  baseline 
code  results  for  a  very  thin  antenna  having  a  radius  change  of  ra 
10  to  1. 


the  wire  of  smaller  radius,  whereas  the  current  computed  by 
NEC  on  the  larger  radius  wire  is  in  excellent  agreement  with 
that  obtained  from  PEC0.  The  charge  distribution  on  the  wire 
with  the  smaller  radius  also  shows  a  slight  difference.  Figs. 
3.5  and  3.6  show  the  computed  current  and  charge  distributions 
which  result  when  the  radius  of  the  larger  wire  is  increased 
to  k  a  ^  *  0.01.  There  is  now  a  significant  difference  in  the 
current  on  the  section  of  the  wire  with  smaller  radius,  but 
on  the  section  with  larger  radius  (which  contains  the  source) 
the  agreement  is  still  excellent.  (Note  that  although  an 
input  admittance  comparison  alone  might  show  the  methods  to  be 
in  excellent  agreement,  comparison  of  the  current  and  charge 
distributions  shows  that  they  are  not.)  Figs.  3.7  and  3.8 
show  current  and  charge  distributions  for  the  same  structure 
calculated  with  NEC  when  a  very  large  number  of  unknowns  (48) 
is  used.  The  differences  between  results  calculated  by  NEC 
and  PEC0  are  even  greater  in  these  figures.  The  use  of  so 
many  unknowns  is,  of  course,  impractical  on  such  a  short 
antenna  under  normal  conditions  but  here  serves  to  illustrate 
a  point.  Upon  close  inspection  of  the  charge  distribution 
calculated  by  NEC  (Fig.  3.8),  one  observes  the  same  type  of 
behavior  in  the  junction  region  as  was  observed  with  our  thin- 
wire  approach  in  Section  II;  namely,  the  charge  seems  to 
attempt  to  satisfy  the  true  edge  behavior  at  the  junction. 

On  the  other  hand,  the  Wu-King  charge  jump  condition  has  also 
been  enforced  and  the  two  conditions  are  incompatible.  One 
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;  .  ,  i  i  hi'iiisi'l  vt'H  as  erroneous  c  u  r  r  e  n  l  s  wIhmu-vc  r  m  a  t  Ii  p  o  i  n  i  ■  ■ 
tor  N  K C  wi  l  e  I  ih  .i  t  I'd  within  a  icw  radii  (  e  .  g  i  i  v  <•  times  the 
radius  of  the  larger  wire)  from  the  junction.  I'nliirtuinli-lv, 
as  reference  to  Figs.  3.3  and  3.4  shows,  the  region  in  which 
effects  of  the  junction  (edge  conditions)  on  the  current  and 
charge  seem  to  come  into  plav  can  be  much  greater  that  this. 
Observation  of  a  great  deal  of  data  has  indicated  to  us  that 
the  actual  distance  from  the  junction  at  which  such  rf I  rets 
become  negligible  seems  to  depend  on  a  number  of  factors  siu  ! 
as  the  wire  lengths,  wire  radii  and  their  ratio.  Broadly 
speaking  ,  however,  at  a  charge  maximum  the  extent  of  the  edgi. 
effects  region  appears  to  be  a  maximum  as  well.  A  <  eminent  i: 
in  order  concerning  the  appearance  of  t1"'  error  onlv  on  the 
thinner  wire  in  these  and  subsequent  figures. 

The*  d  i.  f  f  e  r  e  n  c  e  in  the  current  distrihut  i  o  n  s  o  b  t  a  i  n  e  d  h  v 
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^2  +  1*  •  We  attribute  the  appearance  of  the  "error"  current 
only  on  the  thinner  wire  in  Figs.  3. 5-3. 8  and  in  many  of  the 
subsequent  figures  merely  to  our  choice  of  a  half-wavelength 
length  (  k  L  ^  =  3.05)  for  the  thinner  wire. 

Next,  let  us  consider  two  cases  for  which  the  radius  of 
the  smaller  wire  is  increased  to  ka^  =  0.001.  Figs.  3.9  and 

3.10  show  the  current  and  charge  distributions  on  a  structure 
with  ka^  =  0.00125,  representing  very  small  change  in  radius. 
The  agreement  between  PEC0  and  NEC  is  again  excellent.  Figs. 

3.11  and  3.12  compare  results  when  the  radius  of  the  larger 
wire  is  increased  to  ka^  =  0.01.  For  this  10  to  1  change  in 
radius,  the  current  distribution  calculated  using  NEC  begins 
to  differ  significantly  from  that  of  the  baseline  code  on 

the  segment  with  the  smaller  radius.  (Note  again  that  agree¬ 
ment  is  excellent  on  the  segment  with  the  larger  radius.) 

The  a'-ove  results  are  typical  of  a  large  number  of 
observed  cases.  One  might  expect  that  for  thicker  wires, 
less  radical  changes  in  radii  can  be  made  before  edge  effects 
contaminate  NEC  current  distributions.  Figs.  3.13  and  3.14 
show  current  and  charge  distributions  on  a  rather  thick  wire 
for  a  very  small  change  in  radius:  k  a  ^  =  0.01,  k  a  ^  =  0.0125. 
The  agreement  between  NEC  and  PEC0  is  excellent.  However, 
when  the  radius  of  the  thicker  segment  is  increased  to 
k  a  ^  -  0 . 0  5 ,  one  obtains  the  poorly  agreeing  distributions 
shown  in  Figs.  3.15  and  3.16.  Thus  on  thicker  wires, 
erroneous  currents  appear  with  NEC  for  less  severe  radius  changes. 

Consider  next  the  same  wire  as  in  Fig.  3.11  but  with 
the  radius  of  the  thicker  wire  increased  to  k  a  ^  =  0.05  (i.e.. 
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code  results  for  a  thin  wire  having  a  radius  change  of  ratio 
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Figure  3.10.  Comparison  of  calculated  charge  distribution  (NEC)  with  baseline 
code  results  for  a  thin  wire  having  a  radius  change  of  ratio 
1.25  to  1. 
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Figure  3.12.  Comparison  of  calculated  charge  distribution  (NEC)  with  baseline 
code  results  for  a  thin  wire  having  a  radius  change  of  ratio  10 


Figure  3.14.  Comparison  of  calculated  charge  distribution  (NEC)  with  baseli 
code  results  for  a  fairly  thick  wire  having  a  radius  change  of 
ratio  1.25  to  1. 
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the  radius  change  at  tlie  junction  is  a  ratio  of  5  0  to  1). 

Figs.  3.17  and  3.18  show  the  current  and  charge  distributions 
o  b  t  a i n  e  d  for  t  h i s  case.  As  we  h  a  v  e  seen  before  currents  in 
Fig.  3.17  are  significantly  in  error  ( -r  50%  error)  on  the 
st  ’men t  of  wire  with  the  smaller  radius,  but  again  agree  well 
with  the  baseline  code  results  on  the  segment  of  larger 
radius  (which  contains  the  feed  point).  If  the  values  of  the 
wire  radii  are  simply  reversed,  i.e.  the  driven  segment 
becomes  the  thin  wire  (  k  a  ?  =  0.001)  and  the  un driven  segment 
becomes  the  thick  wire  (ka^  =  0.05),  the  current  and  charge 
distributions  of  Figs.  3.19  and  3.20  result.  In  the  figures, 
one  observes  that  the  current  calculated  with  NEC  on  the 
segment  with  the  larger  radius  differs  significantly  from 
that  predicted  by  the  baseline  code.  On  the  segment  having 
the  smaller  radius  the  real  part  of  the  current  calculated 
by  NEC  agrees  well  with  the  baseline  code  result,  while  the 
imaginary  part  shows  a  slight  shift  toward  the  source  region. 
Thus  it  is  not  true,  as  one  might  be  led  to  suspect  from  the 
previous  results,  that  the  error  tends  to  be  localized  on  the 
thinner  wire.  As  we  have  noted  earlier,  the  error  current 
distribution  probably  depends  on  the  radii  and  lengths  of 
the  two  segments  as  well  as  the  source  location. 

For  all  the  cases  considered  thus  far,  the  junction  has 
been  located  such  that  a  charge  maximum  occurs  at  or  near  the 
junction  region  (discounting  the  region  very  near  the  ju no¬ 
tion  where  an  inf inite  charge  occurs)  .  We  h  a  v  e  purposely 
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Figure  3.17.  Comparison  of  calculated  current  distribution  (NEC)  with  baselin 
code  results  for  a  stepped-radius  antenna  having  a  radius  change 
of  ratio  50  to  1,  where  the  source  is  located  on  the  section  with 
the  larger  radius. 
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Figure  3.19.  Comparison  of  calculated  current  distribution  (NEC)  with  baseline 
code  results  for  a  stepped-radius  antenna  having  a  radius  change 
of  ratio  50  to  1,  where  the  source  is  located  on  the  section 
with  the  smaller  radius. 


limited  ourselves  to  such  cases  to  illustrate  the  influence 
of  the  junction  region  primarily  because  it  has  been  found 
that  an  accurate  treatment  of  the  charge  is  more  critical 
near  a  charge  maximum.  To  illustrate,  consider  Figs.  3.21 
and  3.22,  which  show  the  current  and  charge  distributions  on 
a  s tepped- radius  antenna  structure  for  which  the  junction 
is  located  near  a  charge  minimum  (again,  discounting  the 
region  in  the  neighborhood  of  the  junction).  The  radii  are 
the  same  as  in  Fig.  3.17  (ka^  *  0.001,  ka2  «  0.05)  so  that 
there  is  a  50  to  1  change  in  radius  at  the  junction.  The 
lengths  of  the  various  wire  sections  are  kL^  »  4.7,  kL2  *>  3.2, 
and  kL^  =  0.6.  We  note  that  although  there  is  some  difference 
between  results  calculated  using  NEC  and  those  of  the  baseline 
code  PEC0,  the  differences  are  much  less  significant  than 
those  exhibited  in  Figs.  3.17  through  3.20,  where  the  same 
radii  but  different  wire  lengths  are  employed. 

We  have  also  considered  cases  for  which  the  junction  is 
placed  near  the  feed  point  or  near  the  end  of  the  antenna. 
These  results  have  generally  shown  the  same  trends  as  those 
presented  in  Figs.  3.1  through  3.18  except  that  the  current 
distributions  calculated  by  NEC  appear  to  be  affected  slightly 
more  than  for  a  similar  change  in  radius  occuring  away  from 
ends  or  sources.  The  current  and  charge  distributions  for 
two  sample  cases  are  presented  in  Figs.  3.23  through  3.26. 

For  these  two  cases  the  step  change  in  radius  is  located 
approximately  0.1X  from  the  source.  Observations  of  a  number 
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3.21.  Comparison  of  calculated  current  distribution  (NEC)  with  baseline 
code  results  for  a  stepped-radius  antenna  having  a  radius  change 
of  ratio  50  to  1,  where  the  change  in  radius  occurs  near  a  charge 
minimum. 
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are  3.22.  Comparison  of  calculated  charge  distribution  (NEC)  with  baseline 
code  results  for  a  stepped-radius  antenna  having  a  radius  change 
of  ratio  50  to  1,  where  the  change  in  radius  occurs  near  a  charge 
minimum. 
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Figure  3.23.  Comparison  of  calculated  current  distribution  (NEC)  with  baseline 
code  results  for  a  very  thin  antenna  having  a  radius  change  of 
ratio  10  to  1  approximately  one-tenth  wavelength  away  from  the 
feed. 
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ure  3.24.  Comparison  of  calculated  charge  distribution  (NEC)  with  baseline 
code  results  for  a  very  thin  antenna  having  a  radius  change  o 
ratio  10  to  1  approximately  one-tenth  wavelength  away  from  the 
feed. 
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of  cases  show  that,  if  junctions  are  at  least  0.1A  from 
wire  ends  or  sources,  the  error  in  current  distributions 
computed  by  NEC  is  approximately  the  same  as  when  the  junc¬ 
tion  is  placed  near  a  charge  maximum  away  from  ends  and 
sources  for  similar  radii  wires.  Cases  in  which  the  wire 
with  the  smaller  radius  is  fed  very  near  the  radius  change 
are  of  considerable  practical  interest  and  are  considered  in 
more  detail  later. 

Because  of  the  large  number  of  parameters  on  which 
accuracy  depends,  the  data  presented  here  is  not  easily 
summarized.  Nevertheless,  we  attempt  to  reduce  the  data  to 
a  more  compact  form  from  which  we  derive  a  "rule  of  thumb" 
for  ratios  of  radii  which  NEC  can  handle.  Tc  uo  this  we 
consider  only  current  distributions  and  compare  the  calcula¬ 
tions  of  NEC  and  the  baseline  code  PEC0.  As  an  error 
criterion,  we  choose  the  maximum  magnitude  of  the  difference 
in  current  distributions  calculated  using  NEC  and  PEC0, 
excluding  the  junction  and  feed  regions.  Figs.  3.27  through 
3.29  depict  the  maximum  error  in  the  current  as  a  function 
of  the  radius  change  ratio  a  /a^  where  a^  is  fixed  for  each 
figure.  Separate  curves  are  plotted  for  various  subdomain 
sizes  used  to  model  the  wire  with  NEC.  The  subdomain  sizes 
indicated  in  the  figures  represent  the  width  of  the  smaller 
of  the  two  subdomains  on  each  side  of  the  junction.  Exclud¬ 
ing  curves  for  the  junction  near  the  feed,  all  curves  were 
obtained  for  an  antenna  with  lengths  kL^  =  3.05,  kL^  =  3.2, 
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Maximum  error  in  the  current  computed  by  NEC  for  a  stepped-radius  antenna 
plotted  as  a  function  of  radius  change  ratio  and  subdomain  size  for  a 
reference  wire  with  radius  ka,  =0.0001. 


Figure  3.28.  Maximum  error  in  the  current  computed  by  NEC  for  a  stepped-radius  antenna 
plotted  as  a  function  of  radius  change  ratio  and  subdomain  size  for  a 
reference  wire  with  radius  ka ,  =0.001. 


Figure  3.29.  Maximum  error  in  the  current  computed  by  NEC  for  a  stepped-radius  antenna 
plotted  as  a  function  of  radius  change  ratio  and  subdomain  size  for  a 
reference  wire  with  radius  ka,  =0.01. 


and  kL^  =  0.6,  since  this  results  in  a  r 6  ■  he  maximum  near 
the  junction  (which  gr.itially  makes  the  error  worse).  While 
these  figures  should  not  be  viewed  as  representing  absolute 
error  for  all  possible  wire  configurations,  they  should  con¬ 
stitute  a  rough  upper  bound  to  the  error.  Furthermore, 
general  trends  and  orders  of  magnitude  can  be  deduced  from 
the  figures  with  some  confidence. 

As  the  figures  show,  the  percentage  of  error  (compared 
to  the  baseline  code)  in  current  distributions  calculated 
with  NEC  appears  to  depend  on  three  major  factors:  the  di¬ 
mensions  in  terms  of  wavelengths  of  the  wire  radii,  the  ratio 
of  the  radius  change,  and  the  distance  from  the  junction  to  the 
nearest  match  point  (i.e.,  A/2).  The  distance  of  the  junction 
from  feed  points  and  wire  ends  also  appears  to  be  a  factor 
but  of  less  significance.  We  note  that  in  general  the  error 
decreases  as  the  subdomain  size  increases,  or  equivalently, 
as  the  distance  from  the  junction  to  the  nearest  match  point 
increases,  assuming,  of  course,  that  the  subdomain  widths 
remain  small  enough  so  that  the  current  outside  the  junction 
region  can  be  adequately  modeled.  Our  results  indicate  that 
subdomain  widths  on  the  order  of  0.14\  (kA  ~  0.867)  are  gen¬ 
erally  adequate.  (Note,  however,  that  the  curves  corresponding 
to  this  subdomain  width  in  Figs.  3.27  through  3.29  are  more 
difficult  to  generate  because  interpolation  of  the  data  is 
required  and  thus  the  curves  obtained  may  not  be  as  accuarate 
as  those  obtained  using  smaller  subdomains.  Also  note  that 
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reference  [11]  recommends  the  use  of  subdomains  of  approxi¬ 
mately  0.1A  width  or  less.)  From  the  information  in  these 
figures  one  may  obtain  a  general  idea  of  the  parameter 
ranges  for  which  the  NEC  code  yields  acceptable  results.  If 
a  specific  subdomain  size  is  chosen,  say  kA  -  0.867,  the  data 
in  Figs.  3.27-3.29  may  be  reduced  to  a  more  readily  usable 
form  such  as  Fig.  3.30.  From  this  figure  one  can  quickly  esti¬ 
mate  acceptable  limits  on  the  wire  radii  for  a  fixed  maximum 
allowable  error  in  the  current.  If,  for  example,  the  maximum 
error  is  specified  as  10%,  then  for  wires  with  the  smaller 
radii  ka ^  =  0.01,  0.001,  and  0.0001,  acceptable  limits  on 
the  wires  of  larger  radii  are  ka2  =  0.035,  0.0089,  and  0.0038, 
respectively.  Again,  we  emphasize  that  the  information  sum¬ 
marized  in  this  figure  represents  only  rough  upper  bounds  on 
the  error  of  currents  calculated  by  NEC.  For  many  wire  con¬ 
figurations  and  excitations  the  error  in  the  current  may 
possibly  be  larger,  but  is  generally  expected  to  be  less 
significant. 

One  may  derive  a  rough "rule  of  thumb"  formula  for  ac¬ 
ceptable  ratios  of  radii  (i.e.,  error  less  than  5%)  by 
requiring  that  all  combinations  of  radii  fall  below  the 
straight  line  drawn  in  Figure  3.30.  Thus,  if  the  closest 
match  point  is  about  0 .  14X  from  the  junction,  radii  ratios 
which  satisfy  the  condition 

*'7^  -  ~  ~h  Un(kal)  +  3]  (3.2) 


with  h_  ai  can  &e  expected  to  yield  reasonably  accurate 
current  distributions.  It  is  clear  from  the  figure,  that 


for  very  thin  wires  (2,n 


>  10),  the  above  estimate  is 


probably  too  conservative.  The  equation  also  appears  to  be 
too  conservative  for  thicker  wires  ( &n,  —  -  <  3),  but  failure 
of  the  thin-wire  assumptions  more  than  likely  precludes  the 
use  of  NEC  in  this  range  for  most  applications. 

The  results  obtained  to  this  point  apply  to  antennas 
only.  We  next  consider  very  briefly  two  stepped-radius 
scat terers .  Both  structures  are  approximately  one  wavelength 
long  and  have  a  step  change  in  radius  midway  between  the  two 
ends.  Fig.  3.31  illustrates  the  current  distribution  excited 
by  a  normally  incident  plane  wave  when  ka^  =  0.001  and 
ka2  =  0.005.  For  the  same  structure,  Fig.  3.32  shows  the  current 
distribution  excited  by  a  plane  wave  which  is  incident  at  an 
angle  of  45°  with  respect  to  the  positive  l  axis  (see  Fig. 
2.1).  The  phase  reference  in  each  of  the  figures  is  located 
at  s  =  0.  Results  obtained  with  the  NEC  code  used  a  sub¬ 


domain  width  of  kA  =  0.53.  From  Fig.  3.28,  we  expect,  for 
subdomains  of  this  width  and  a  radius  change  ratio  of  5  to  1, 
a  maximum  error  in  the  current  calculated  by  NEC  of  about  8% 
or  less.  This  is  approximately  the  maximum  error  present  in 
Figs.  3.31  and  3.32.  Figs  3.33  and  3.34  show  the  current 
distributions  excited  by  plane  waves  with  normal  and  45° 
incidence,  respectively,  for  a  structure  with  radii 
kaj  =  0.001  and  ka^  =  0.05.  Again,  from  Fig.  3.28  for  a  radius 
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Figure  3.34.  Comparison  of  calculated  current  distribution  (NEC)  with  base 
code  results  for  a  stepped-radius  scatterer  excited  by  an  obi 
incident  plane  wave  and  having  a  radius  change  of  ratio  50  to 


change  ratio  of  50  to  1,  we  expect  to  find  a  maximum  error 
in  the  current  of  about  45%.  This  predicted  error  corre¬ 
sponds  well  with  the  actual  error  in  the  currents  computed 
by  NEC  in  Figs.  3.33  and  3.34.  In  general,  currents  calcula¬ 
ted  with  NEC  for  scattering  structures  appear  to  exhibit  the 
same  types  and  orders  of  error  as  for  driven  structures. 

This  should  be  expected,  since  excitation  of  a  scatterer  can 
be  thought  of  as  an  antenna  problem  with  a  distributed  source 
Furthermore,  derivation  of  the  charge  jump  condition  does 
not  depend  on  the  source,  except  that  point  sources  located 
directly  at  the  junction  are  excluded. 

We  next  study  an  antenna  driven  either  in  the  annular 
junction  region  (see  Fig.  2.1)  or  on  the  '..ire  of  smaller 
radius  at  a  point  very  near  the  junction.  Such  a  configura¬ 
tion  approximately  models,  for  example,  a  whip  antenna 
driven  against  a  kingpost.  The  baseline  code  PEC0  is  capable 
of  treating  stepped-radius  structures  which  are  driven  a  n  v  - 
where  along  their  length,  including  within  the  annulus  region 
The  NEC  code,  on  the  other  hand,  is  not  presently  capable 
of  treating  structures  which  are  fed  precisely  at  a  step 
change  in  rad ius --becaus e  it  assumes  that  delta  gap  sources 
are  located  at  the  center  of  subdomains,  while  all  changes 
in  radius  occur  at  subdomain  boundaries.  NEC  does  offer  an 
alternative  point  source  model  called  the  "current-slope- 
discontinuity  voltage  source"  [11],  which  is  effectivelv 
located  at  the  boundaries  between  subdomains.  According  to 
[11],  however,  no  change  is  radius  is  permitted  at  the  point 
where  this  source  model  is  used. 
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We  first  attempt  to  use  NEC  to  simulate  a  s tepped-radius 


antenna  fed  at  the  junction.  For  a  delta  gap  source  model, 
this  might  be  accomplished  by  replacing  the  original  junction 
with  a  short  segment  of  wire  of  intermediate  radius.  Figs. 
3.35  and  3.36  depict  the  current  and  charge  distributions, 
respectively,  on  a  1.02A  antenna  which  has  a  step  radius 
change  in  the  center  of  the  antenna  and  is  fed  in  the  annulus 
region.  The  wire  radii  are  ka1  -  0.001  and  ka2  =  0.005. 

NEC  is  used  to  simulate  the  structure  by  using  a  short  segment 
(kA  *  0.5)  of  radius  ka  »  0.003  in  place  of  the  junction. 

We  note  that  currents  calculated  with  NEC  are  significantly 
in  error.  We  attribute  the  error  exhibited  by  NEC  in  this 
case  partially  to  the  differences  in  the  two  assumed  models 
and  partially  due  to  the  increased  sensitivity  of  the  cur¬ 
rent  distribution  to  such  differences  because  of  the  near 
half-wavelength  wire  elements  chosen. 

If  the  feed  point  is  next  chosen  on  the  wire  with  the 
smaller  radius  and  a  short  distance  from  the  junction,  NEC 
can  apparently  model  the  structure  quite  accurately.  Figs. 
3.37  and  3.38  show  current  and  charge  distributions  obtained 
when  the  feed  is  located  0.0424A  away  from  the  junction.  The 
agreement  with  the  baseline  code  for  this  case  is  excellent. 
To  move  the  feed  point  closer  to  the  junction  requires  a 
reduction  in  subdomain  widths  and  hence,  since  the  charge  is 
then  defined  closer  to  the  junction,  we  might  expect 
the  errors  in  the  currents  calculated  by  NEC  to  be  somewhat 
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Figure  3.35.  Current  distribution  on  a  stepped-radius  antenna  which  is  fed  in 
the  annulus  region.  For  the  NEC  code  the  junction  is  simulated 
with  a  short  wire  of  intermediate  radius  ka  =  0.003. 
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Figure  3.36.  Charge  distribution  on  a  stepped  radius  antenna  which  is  fed  in 
the  annulus  region.  For  the  NEC  code  the  junction  is  simulated 
with  a  short  wire  of  intermediate  radius  ka  =  0.003. 


Comparison  of  calculated  current  distribution  (NEC)  with  bas 
code  results  for  a  stepped-radius  antenna  fed  on  the  wire  of 
smaller  radius  and  near  the  junction. 
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results  for  a  stepped-radius  antenna  fed  on  the 
er  radius  and  near  the  junction. 


greater.  Indeed,  Figs.  3.39  and  3.40  confirm  this  conjecture 
for  a  feed  located  0.0127A  away  from  the  junction.  Referring 
to  Fig.  3.28  again  for  k  A  =  0.16  and  a  5  to  1  radius  change 
ratio,  we  might  expect  a  maximum  error  of  about  14%  in 
the  current.  Fig.  3.39  indicates  that  the  currents  calcu¬ 
lated  by  NEC  agree  with  the  baseline  code  results  to  about 
10%.  A  comparison  of  the  current  and  charge  distributions 
in  Figs.  3.37  through  3.40  with  Fig.  3.35  and  3.36  also 
reveals  that  as  the  feed  point  moves  closer  to  the  junction, 
the  distributions  approach  those  obtained  with  the  baseline 
code  for  the  case  where  the  feed  is  placed  directly  in  the 
annulus  region. 

Before  concluding  this  section  we  note  one  further 
interesting  fact  regarding  NEC's  capability  to  model  a 
source  in  the  junction.  We  refer  again  to  the  case  in  which 
the  source  is  placed  in  the  annulus  region  and  note  that,  if 
a  radius  change  were  allowed  across  the  current -slo pe¬ 
dis  continuity  voltage  source  in  NEC,  this  case  could  hi' 
modeled  bv  NEC.  If  we  merely  ignore  the  warning  in  111]  to 

the  radii  on  each  side  of  the  current-slope-discontinuity 
source  model  the  same,  the  current  and  charge  distributions 
of  Figs.  3.41  to  3.44  are  obtained.  While  some  magnitude 
errors  are  apparent  in  the  current  distributions  as  calculated 
bv  NEC,  the  errors  are  within  those  predicted  in  Fig.  3.28 
for  junctions  of  the  radii  used  in  the  figures.  It  appears, 
therefore,  that  the  potential  exists  in  NEC,  possibly  with 
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code  results  for  a  stepped-radius  antenna  which  is  fed  very 
the  junction  region  on  the  wire  of  smaller  radius. 
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code  results  for  a  stepped-radius  antenna  fed  in  the  annulus 
region.  The  current-slope-discontinuity  voltage  source  has 
used  (improperly)  to  obtain  the  results  with  NEC. 
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w , 


pome  modification  and/or  justification  of  the  present  treat¬ 


ment  of  the  current-slope-discontinuity  voltage  source,  to 
accurately  model  structures  which  are  fed  at  junctions. 
Alternatively,  it  may  be  possible  to  modify  the  Wu-King  con¬ 
dition  so  as  to  incorporate  a  voltage  source  located  at  the 
junction.  In  the  Appendix,  it  is  argued  that  the  Wu-King 
junction  condition  can  be  interpreted  as  a  condition  on 
continuity  of  scalar  potential  across  the  junction.  If  a 
voltage  source  is  present  at  the  junction,  the  condition 
should  be  modified  so  that  the  scalar  potential  is  discon¬ 
tinuous  by  the  voltage  V  of  the  source  at  the  junction.  Thus 
from  Eqs.  (A.3)and  (A. 4)  of  the  Appendix  we  obtain 

27rev  =  q 2  t£n(^-)  -  y]  -  qx  Un(j^-)  -  y],  (3.3) 

where  the  positive  reference  potential  of  the  source  is 
taken  to  be  on  the  wire  of  radius  a^.  By  the  continuity 
equation,  the  corresponding  current  slopes  at  the  junction 
are 


-jwqj 


9 


-jwq2 
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from  which  a  current-slope-discont inuity  voltage  source  is 
obtained  satisfying 


It  is  interesting  to  compare  (3.4)  with  the  present  current- 
slope-discontinuity  voltage  source  used  in  NEC  when  the  radi 
on  each  side  of  the  junction  are  the  same,  i.e.,  =  a 
In  NEC,  this  source  condition  is  equivalent  to 


V 


jn 

27rk 


C  £  n  (— ) 
a 


i] 


(3.5) 


where  d  is  the  half-length  of  the  subdomains  between  which 
the  source  is  located,  whereas  (3.4)  reduces  to 


V 


y] 


(3.6) 


It  is  seen  that  only  when  we  choose 


kd  =  e 


1.526 


do  the  above  two  equations  agree.  It  is  claimed  in  [111, 
however,  that  good  results  are  obtained  with  Eq .  (3.5) 
reqardless  of  the  value  of  d.  On  the  other  hand,  (3.5)  does 
not  show  the  logarithmic  dependence  on  frequency  evident  in 
(3.6)  which  comes  from  (3,3)  and  has  already  been  shown  to 
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be  correct  in  the  form  of  Eq .  (2.12).  We  have  not  managed 

to  resolve  these  Inconsistencies  during  this  study--indeed , 
we  have  derived  still  other  expressions  differing  from  both 
(3.5)  and  (3.6)  only  in  the  form  of  the  logarithmic  terms. 
This  should  not  be  too  surprising,  however,  since  it  is  well 
known  that  in  making  approximations  in  linear  antenna  theory, 
constants  associated  with  the  terms  involving  logarithms 
of  the  radius  cannot  be  uniquely  fixed.  We  leave  unresolved 
at  this  time  the  problem  of  combining  the  two  current-slope- 
discontinuity  conditions— one  for  the  s tepped-radius ,  and 
one  for  sources-- into  a  single  comprehensive  condition. 
However,  we  suggest  that  this  may  be  a  fruitful  area  for 
further  research. 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  report  we  have  considered  methods  for  tieating 
junctions  of  collinear  wires  with  dissimilar  radii.  A  nu¬ 
merical  technique  has  been  presented  for  accurately  computing 
current  and  charge  distributions  on  a  stepped-radius  antenna. 
The  method,  in  fact,  yields  the  current  and  charge  distribu¬ 
tions  in  the  annular  junction  region,  permitting  one  to  study 
their  behavior  there  in  detail.  The  technique  employs  a 
body  of  revolution  model  for  the  stepped-radius  wire  and 
hence  treats  only  junctions  of  collinear  wires  .  An  ancillary 
result  of  this  study  has  been  the  development  and  validation 
of  a  thin-wire  theory  approach  which  can  also  be  used  to 
accurately  calculate  current  and  charge  distributions  away 
from  the  junction  region.  The  thin-wire  approach  retains 
the  advantages  of  the  thin-wire  approximations,  permits  the 
use  of  the  very  simple  pulse  basis  set,  and  does  not  require 
any  condition  to  be  enforced  on  the  charge  at  the  junction. 
Furthermore,  since  it  has  been  demonstrated  that  no  condition 
is  required  on  the  charge  at  junctions  it  is  conjectured  that 
the  thin-wire  approach  is  capable  of  accurately  treating 
junctions  of  non-co 1 1  inear  wires  of  dissimilar  radii  and  mul¬ 
tiple  wire  junctions. 


1 1  7 


The  Nw.or  u'.t  !  V  l  tH' t  v  om  .-igne  t  i  v  -;  (\>de  t\Ff'  u  ;  .•  ->  *■  i  ,  i 
functions  which  requite  .1  c  etui  1  (  1  >•  n  »<n  1  It .  .it  .it  ,-  11 

boundaries.  For  segments  of  the  same  radii,  this  ooudil  ion 
is  merely  that  the  charge  be  continuous  from  one  segment  to 
the  next.  When  adjacent  segments  have  different  radii, 
however,  this  condition  is  no  longer  valid.  To  test  the 
validity  of  various  proposed  charge  "jump  conditions"  at  a 

radius  discontinuity  both  the  body  of  revolution  and  the  thin- 
wire  codes  have  been  used  to  calculate  values  of  the  extrapo - 
lated  charge  on  both  sides  of  a  radius  step,  from  which  a 
numerical  charge  jump  condition  can  be  calculated.  All  our 
numerical  results  have  verified  that  the  analytical  charge 
jump  condition  of  Wu  and  King  [7]  is  indeed  satisfied  at  junc¬ 
tions  of  collinear  wires.  At  the  present  time,  no  accurate 
model  exists,  nor  does  one  appear  feasible  to  develop,  to 
verify  the  Wu-King  condition  for  non-collinear  or  multi-wire 
junctions  of  unequal  radii.  However,  we  conjecture  that  the 
condition  is  also  correct  for  these  configurations.  We  fur¬ 
ther  suggest  that  the  thin-wire  approach  could  be  extended  to 
these  configurations  and  should  yield  charge  distributions  sat¬ 
isfying  the  Wu-King  condition.  This  is  because  the  thin-wire 
approach  requires  (approximately)  continuity  of  scalar  poten¬ 
tial  at  wire  junctions  and  this  condition  is  also  (approximately) 
enforced  by  the  Wu-King  condition.  The  authors  suggest  that  such 
a  test  of  the  thin-wire  theory  would  be  useful  as  a  check  on  both 
the  validity  of  the  thin-wire  approach  and  on  the  Wu-King 
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junction  condition  for  multi-wire  junctions.  Although  the 
Wu-King  charge  jump  condition  is, strictly  speaking,  valid 
only  at  least  several  radii  from  the  junction,  the  NEC 
computer  code  [11]  explicitly  enforces  the  condition  at  all 
wire  junctions.  The  accuracy  of  the  code  has  been  studied 
for  various  locations  of  the  step  with  respect  to  feed  points 
and  free  wire  ends  and  with  respect  to  its  location  at  charge 
minima  or  maxima.  The  results  obtained  by  comparing  the  cur¬ 
rent  and  charge  distributions  calculated  with  NEC  co  those 
computed  with  the  body  of  revolution  code  generally  indicate 
satisfactory  performance  of  NEC  provided  certain  limitations 
are  observed.  Three  primary  factors  have  been  found  to  influ¬ 
ence  the  accuracy  of  results  obtained  with  NEC:  the  radius  of 

the  smaller  wire  in  terms  of  wavelengths,  the  ratio  of  the 
change  in  radius,  and  the  distance  from  the  junction  to  the 
nearest  match  point.  Information  has  been  presented  in 
Section  III  which  should  enable  users  to  determine  appropriate 
limitations  to  impose  on  subdomain  sizes  and  radius  change 
ratios.  Furthermore,  an  empirical  formula  is  given  which 
relates  the  maximum  allowable  radius  change  to  the  value 
of  the  smaller  radius  in  terms  of  wavelengths  whenever  the 
subdomain  size  is  restricted  to  about  1/7A.  It  may  be  desir¬ 
able  to  incorporate  this  condition  into  the  NEC  manual  and/or 
to  have  NEC  print  warning  messages  to  the  user  whenever 
the  condition  is  not  met. 
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It  is  recommended  that  some  consideration  be  given  to 
attempting  to  combine  the  Wu-King  type  junction  condition 
with  a  s 1 op e- d i s c on t  i  nu  i  t y  source  condition.  Both  condi¬ 
tions  on  the  charge  discontinuity  arise  from  considering 
the  difference  in  potential  between  adjacent  segments.  Some 
of  the  cases  considered  in  this  study  indicate  that  combining 
the  two  conditions  might  permit  NEC  to  handle  a  source  placed 
directly  at  a  stepped-radius  junction. 
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APPENDIX 


In  this  appendix  we  derive  the  Wu-King  jump  condition 
[7]  for  the  charge  near  a  s  tepped-radius  junction.  The 
derivation  given  here,  while  less  rigorous  than  that  given 
in  [7],  is  simple  and  yet  convincing.  It  is  shown  that  for 
thin  wires  the  condition  is  equivalent  to  requiring  that 
the  potential  be  continuous  across  the  junction. 

Consider  the  stepped-radius  antenna.  Fig.  2.1,  where 
it  is  assumed  that  the  wire  radii  a^  and  a ^  are  both  elec¬ 
trically  small  (i.e.,  ka^,  ka^  <<  1).  Now  consider  the 

evaluation  of  the  scalar  potential  at  a  point  Aj  units 
from  the  junction  on  the  axis  of  the  wire  of  radius  a^. 

Only  the  circumferentially  uniform  component  of  the  charge 
contributes  to  the  potential  on  the  axis,  sc  that  ignoring 
the  small  amount  of  charge  on  the  an nulls  region  connecting 
the  two  wires  of  different  radii  at  the  junction,  we  have 
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where 


2  7  7 

Ri  -  <*i  -  »■)  +  =;  • 

2  7  7 

R2  =  (Sj  -  s’)2  +  a 2  , 

and  where  s  ^  +  L  ^  +  A  ^  ,  and  is  the  scalar  potential 

evaluated  at  s^.  We  assume  that  A^  is  chosen  such  that 

max  ( ka  ^ ,  ka2)  <<  kA^  <<  1.  (A. 2) 

Thus,  while  our  observation  point  is  located  several  radii 
from  the  junction,  it  is  still  near  the  junction  electrically. 
Such  a  point  can  always  be  found  provided  the  wire  radii  are 
sufficiently  small.  Obviously,  most  of  the  contribution  to 
the  integrals  in  (A.l)  comes  from  the  neighborhood  of  the 
point  s'  =  Sj  if  a  ^  is  small.  This  permits  three  approxi¬ 
mations:  (1)  Since  q  (  s  '  )  is  relatively  slowly  varying  in  the 
neighborhood  of  s'  =  s^,  we  approximate  it  as  a  constant 
(recall  that  we  are  sufficiently  far  removed  from  the  junction 
to  be  out  of  the  region  where  edge  effects  dominate);  thus 
q ( s )  -  q(s^)  =  q  ^ .  (2)  For  a^  sufficiently  small,  the 
first  integral  in  (A.l)  can  be  neglected  compared  to  the 
second.  (3)  Finally,  since  most  of  the  contribution  to  the 
second  integral  comes  from  the  neighborhood  of  s'  =  s^,  the 
limits  of  integration  are  unimportant  and  may  be  extended  to 
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infinity  in  both  directions.  While  these  approximations  are 
somewhat  gross,  a  more  rigorous  analysis  leads  to  the  same 
result  as  that  below  at  considerable  expense  in  simplicity. 
With  the  above  approximations,  (A.l)  becomes 
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ds  ' 


4j  e 


(ka1> 


=  Hi 
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-  Y) 


(A. 3) 


where  y  =  0.5772  is  Euler's  constant  and  the  small  agument 
approximation  of  the  Hankel  function,  H^2^(x)  =  -  -^-(in^  +  y) , 
has  been  used. 

Now  evaluate  the  potential  at  a  point  s  =  a  distance 
from  the  junction  on  the  wire  of  radius  a^.  Assume 
also  satisfies  (A. 2).  Then  by  similar  reasoning  as  before, 
the  potential  $(32)  =  $2  is  approximately 


$  2 

2 


2ne 


( in 


ka , 


Y> 


(A. 4) 


Since  the  potential  is  a  smoothly  varying  quantity  and  since 
both  observation  points  are  at  small  (electrical)  distances 
from  the  junction,  the  potentials  in  (A. 3)  and  (A. 4)  are 
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approximately  equal.  Equating  them  yields  the  Wu-King 
condition  [ 7 ]  , 


qi  (£n  kT 


-  Y) 


q2  Un  kT 


Y)  • 


(A.  5) 


Thus  we  see  that  the  Wu-King  condition  may  be  interpreted  as 
arising  from  the  continuity  of  scalar  potential  at  the 
junction.  We  hasten  to  point  out,  however,  that  the  latter 
condition  is  always  true,  even  when  the  wires  become  thick  so 
that  the  Wu-King  condition  no  longer  holds.  Furthermore,  the 
Wu-King  condition  merely  relates  the  behavior  of  the  charge 
near  the  junction,  (recall  Eq .  (A. 2))  but  ignores  the  de¬ 

tailed  behavior  of  the  charge  in  the  junction  region  itself. 

If  is  also  interesting  to  note  that  the  Wu-King  con¬ 
dition  is  also  easily  derivable  from  the  discussion  of 
Schelkunoff  and  Friis  [18],  who  determine  the  current  on 
either  side  of  s tepped-radius  junction.  They  also  determine 
the  behavior  of  the  current  at  the  junction  by  enforcing 
continuity  of  current  and  scalar  potential  there. 
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